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I. INTRODUCTION 



The subject of special functions has a long and varied history. At the core of special 
functions lie the Gaussian hypergeometric function and its confluent forms, the confluent 
hypergeometric functions. The slight modification of the confluent hypergeometric functions 
are also known as Whittaker functions. 

The confluent hypergeometric functions include as special cases Bessel functions, parabolic 
cylinder functions. Coulomb wave functions, and incomplete gamma functions. Numerous 
properties of confluent hypergeometric functions flow directly from a knowledge of the Gaus- 
sian hypergeometric function, and a basic understanding of the Gaussian and the confluent 
hypergeometric functions is sufficient for the derivation of many characteristics of all the 
above named functions. A natural generalization of the Gaussian hypergeometric function is 
the generalized hypergeometric function, which in turn is generalized by Meijer's G-function. 
The theory of the generalized and the G-function is fundamental in the applications, since 
they contain as special cases all commonly used functions of mathematical analysis. Further, 
these functions are the building blocks for many other functions which do not belong to the 
hypergeometric family. 

Recently, taking our lead from the exact solution of a super symmetric differential 
equation^ describing quark-antiquark potential models and their relation to the multi- 
quark systems, namely quark-diquark system and associated potential models (diquark is a 
pair of quarks formed as a point-like particle at one end of a baryon bag) in the description 
of relating the baryon masses to the meson masses, we are led into a discovery of a new 
kind of special function theory that generalizes the confluent hypergeometric series that 
we call Grand Confluent Hypergeometric functions, and show generation of a class of exact 
solutions that may be used in many areas of physics and mathematics and most importantly 
for our own purposes, the supersymmetric particle theories. 

The supersymetric theory we developed earlier--'^ based on supergroups of the type 
SU{m/n) relates baryon masses and meson masses to each other with noteworthy accu- 
racy, but also predicts multiquark states that have been observed experimentally, (such as 
ao and /o mesons made up of diquark-antidiquark pairs) with significant accuracy. In the 
process we have also developed a algebra of color- based on the units of highest division 
algebras, that of octonions that we call "split octonion algebra" that produces the correct 
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color degrees of freedom for quark systems. Split octonion algebra based on the complex- 
ification of the octonionic units is a closed Heisenberg algebra whose properties have been 
described in some of our earlier papers^*^. 

The notion of the solution having to give rise to an extra degree of a quantum number 
we designate as "i*'^ kind of hidden radial quatum number" that will be expressed below, 
makes the theory especially applicable to all supersymmetric theories in nature. Indeed, it 
has an infinite number of eigenvalues, and each one has infinite eigennumbers. For example 
a hydrogen like atom wave function, only has one eigenvalue and has infinite eigennumbers. 
As we see in Regge trajectory plot of angular momentum vs. square of mass (J vs. m^), 
there are many linearly increasing lines with same slopes including bunch of eigenvalues 
corresponding to fermions and bosons. It is not clear what the meaning of many eigenvalues 
are. With our analysis, it will be explained below: Look at fl28ap - fl29|) . The formalism we 
develop in this paper is very general in scope, as it unifies all the special functions mentioned 
above in mathematical physics. It has a remarkable range of applications in mathematics 
and other sciences, as well as large range of applicability in all sorts of engineering problems. 

There is the spin free hamiltonian involving only scalar potential for the q — q system:-!^ 

= 4[(m + hrf + + ^^^] (1) 

where -P^^ = — ^ — f m= mass, b= real positive, and /= angular momentum quantum 
number. If the mass m is negligible because of its extremely small value. ([T]) turns to be 

i72^4[i6V + p2 + i(i±l)] (2) 
When wave function ijjir) = e~t^V'F(r) acts on both sides of ([2]), we have 

Let z = ^br"^ then, ([3]) becomes 

d^F OF 
^^ + (7-^)^-«oF(^) = (4) 

where 7 = / + | and «o = ~ + I) © equivalent to confluent 

hypergeometric equation. Also, ao must be zero or a negative integer. If it doesn't, poly- 
nomial F{z) will blows up as z and n goes to infinity. Physically, we are not interested 
in this kind of situation. The function F{z) must to be finite. Let's set = 1 — iiq 

3 



where no(called primary radial quantum number) = 1,2,3, Then, its eigenvalue is 
ii^^ = 46 [2|ao| + + §)] where |ao| = uq — 1. Also, the normalized wave function by using 
orthogonal relation is 



(26) 


'+§( 


' + !)!(/ + 


1)! 


(2/ + 2)! 


ao 


!(/+ 


ao 





Now let's define normalized wave function including the small mass m. When wave function 
ilj{r) = e"*'-'''^") r^y{r)Yl^* {9,(f)) acts on both sides of ([T]), it becomes 



+ i-hr'^ - 2mr + 2il + 1)) ^ + (Qr - 2mil + 1)) y = 



(6) 



where 



(7) 



By using the function y{r) as Frobinous series in ([6]), we have two indicial roots which are 
Ai = and A2 = —21 — 1. Recurrence formula for all n is 



Kji — 



Bn 



Kn-1 

A, 
B. 



Cn-1 

2m(n+\+l+l) 
" ~ (n+A+l)(n+A+2(/+l)) 

-Q+b{n+\-l) 



(n+A+l)(n+A+2(i+l)) 

Let's investigate function y{r) as n and r go to infinity. As n ^ 1, ([H]) is 

2m b/n 

hm Kn = + -r-^ 

n>i n hm Kn-i 



(8) 



(9) 



We see the first term of r.h.s. in is negligible, since mass m is extremely small and n is 
too large, respectively. Then, (|9]) is approximately equal to 



b/n 



n>i hm Kn-l 

r).>l 



(10) 
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We can classify C„ to its even and odd parts from ( ITOl) by using Kn = %^ and Kn-i = J^" ■ 

Cq Ci = mCo 

C2 = bCo C3 = ImCo 

C4 = Y^b'^Co C5 = 2^mCo 

Ce = ^6^00 C7 = 4mCo (n) 



From (ITT!) suggesting Co = 1 , the function y{r) approximately is 

lim'ufr)^'^- — T^i-br^] +mr'S~^ — i -br^ ] >(l + mr)e2'"' (12) 

It is unacceptable that wave function 'il){r, 6, 0) is divergent as r goes to infinity from the 
quantum mechanical point of view. As r is extremely large value, the big polynomial of 
degree n will take a dominant position. Substitute (fT2l) into the wave function which gives 
ilj{r,6,(f)) = Ne~^^^^~^ r^y{r)Yl^* {d,(f)) ; N is normalized constant. 

limV(r,^,0) > A^r'e~3(^+¥)'(l + mr)e^''"V;'"'(^,0) 

n>l 

^ A^(l + mr)r'e3^"'y;"*(^,0) (13) 

Even if the mass m is extremely small, the wave function ilj{r, 6, 0) will blows up as r— > 00. 
Then it yields 

lim ^(r, e, 0) > lim N{1 + mrye^^''^Yl^* {6, 0) ^ 00 (14) 

r— >oo 

All wave functions must to go to zero as r goes to infinity from a quantum mechanical per- 
spective. The first and second term of y{r) must also be terminated to become a polynomial 
of degree in this case. As we see in (|T2l) . the first term indicates even term of C„, and the 
second term of it has odd term of if y{r) is not terminated at certain value of n, it will 
blow up as r goes to infinity. This is the reason why we argue that y{r) must be terminated 
to become a polynomial of degree n. Now, let's try to define the first kind of independent 
solution as Ai = 0: 



As Ai = 0, 

, , 2{n + l + l)m . , 

= (n + l)(n + 2(/ + l)) ^'''^ 

-Q + bjn - 1) 
= (n + l)(n + 2(/ + l)) 

We define Bij^^^i refering to BiBjE^Bi. We can classify C„ to its even and odd parts up to 

tlie first order of small mass m from ([8]), fllSal) and fllSbl) . . 



Co 




Ci = 


CoAo 




C2 = 


CqBi 


C3 = 


CoiAoB^ + A^Bi) 






CqBi-^ 


c,= 


Co{AqB2^4 + A2B14 + AiBi^-i) 




c,= 


Co-Bl,3,5 


C7 = 


Co(^0-B2,4,6 + ^2-Bl,4,6 + ^4-Bl,3,6 + ^6-Bl,3,5) 








(16) 


Cs = 


Co-Bl,3,5,7 




Co (^0-82,4,6,8 + ^2-Bl,4,6,8 + ^4-Bl,3,6,8 + ^6-Bl,3,5,8 


+ ^8-81,3,5,7) 


Cio = 


= Co-Bl,3,5,7,9 




= 6*0(^0-82,4,6,8,10 + ^2-81,4,6,8,10 + ^4-81,3,6,8,10 





+^6-81,3,5,8,10 + ^8-81,3,5,7,10 + ^10-81,3,5,7,9) 

As we describe y{r) as a power series by using (1161) . 

00 00 

y{r) = y{r)domin. + y{r)small = ^^an^^" + ^ Csn+l'^^"^^ 

n=0 n=0 

1 + 5^ n ^2fe+ir'" [ + CoMor + iAoB2 + A2Biy 

n=l k=0 ^ ^ 

00 p n—1 n— 1 ^ rt— 1 J — 1 \ n—l -i 

+ ^0 n ^2p+2 + Yi (^2, n ^2p+2 n ^2^+1 ) + ^2^ n ^^p+i [ (17) 

n=2 L p=0 j=l ^ p=j A:=0 ^ p=0 ^ 

We clioose one of -82^+1 to be zero, where = 0, 1, 2, 3, ■ ■ ■ in order to make a polynomial 
of degree n in y{r)domin. of (fT6l) . In other words, we might choose Q = 2b{no — 1) where 
uq = 1,2,3,--- in fllSbp . Then -82/0+1 will be zero at certain value of k. However, since 



yij) small in (HTI) : look at the odd term of C„ in (ITBi) . it has combinations of -B2fc-f2and B2k+i 
for every each of the odd term of C„. It means that one of each of -B2A:+2and -82^+1 terms 
must be zero at same time. In other words, Q = 2h{nQ — l)and Q = h{2ni — 1) where 
rii = 0, 1,2,3,--- must be satisfied in this series spontaneously. We call no as primary 
radial quantum number and rii as the first kind of hidden radial quantum number. Then 
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y{r) will be a polynomial degree of n. y{r) consists of two terms which are y{r)domin. and 
yij') small in f|T7|) . The dominant wave function y{r)domin. which does not include small mass 
m must be terminated to become a polynomial of degree n in this case. But the y{r)smaii is 
extremely small wave function because it includes A„|a=o = (n+ixti.+2\'l+i)) ^hich has small 
mass m. Also (I7j) is equivalent to 26(no — 1) and h{2ni — 1) at same time. And in this article 
Pochhammer symbol (x)^ is used to represent the rising factorial: = ^^y{x)^ ■ 
For simplicity, plugging Co = ^^^n^iS'' in ( fT7|) . we obtain 

3 1 

y{r) = Q>V|ao|,|ai|(|ao| = ^0 - 1, |ai| = ni - l,-f = I + -] z = -br'^ 

(7; z) only if |ao| < l^il (18) 

|ao| 



where, 



r(|ao| +7) ^ i-\ao\)n ^n f.Q^ 



yrM r(|ao|+7) HaoDn ^n ' v^" + |(7 " |))r(n + |)r(7 + n - - \a,\), ^^ 

11^1 l^'^^ r(7) ^ n!(7)„ ^ ^ r(A; + n+|)r(A; + n + 7 + i) 



27TiB{\ai\ + 1, 

And, 



mi 1^1 ^ r(|ai| + i)r(i/2) 

"^('"^' + ^'2^= r(K| + 3/2) (2°^) 

«;, = ^^(l-p^) (20b) 
in the above, T(s, u) is the operator which acts on 

Zu(l — t) /-, 9\ 

r{s,t,p,u)= ds s-'/\l + sp^'^'^+^Z^^ dtf'-''/^ dpidu^-— (21) 



^0 j-i 



We see in f llQap . it is the first kind of confluent hypergeometric polynomial of degree |ao| 
surprisingly. (ITSll denoted as QW|o(,|jai| ^|ao| = rio — 1, |ai| = rii — 1, 7 = / + |; z = |&r^j 
is called as the first kind of Grand confluent hypergeometric polynomial of degree |ao| and 
I ail with the first order m. 

Also, we obtain two eigenvalues which are Eq = 46|/+2r2o — | and Ef = 46|/+2ni + l|'. 
The former is called as primary radial eigenvalue. And the latter is called as the first kind 
of hidden radial eigenvalue. As we let the small mass m goes to zero in ( fT8|) . its solution is 



same as the 1'^* kind of confluent hypergeometric polynomial. From ( IT8|l . the wave function 
of it is 

By using orthogonal relation, normalized constant N is 
N-- 



M nrn ^ (-l)l"«l2^-i f 2r(|«o| + 7 - |)r(|ao| + 7 + |) 

r(|ao| + l)r(|ao| +7) 



67 ■ ■ - ^ r(7-| 

|ao| |ai|-n 



EE 



ri=0 fc=0 

where, 



(n + l(7-i))r(|ao 


+ 7)r(7 + 




Oio\)n{n - 




r(7)r(A; + r2- 




+i)(7)n( 







(23) 



|ao| 


= ^0 


- 1 




= ni 





only if |ao| < |«i| (24) 

= Til — 1 I 

([22]) is called the I'** kind of the normalized Grand confluent hypergeometric wave function 
of degree |ao| and |ai|. 

As the small mass m goes to zero in ( l22ll and ( l23il . it turns out 



^^Q-fv , ,H^jno,m,i,m y 2T-ir(|ao| + l)r(|ao| +7) J I K 'rj K J 

(!25|) is exactly equivalent to ([5]). We see in the new special function, there are two eigen- 
numbers uq and rii. The flrst eigennumber called primary radial quantum number appears 
in zeroth order of m term,?/(r)rfomm. and flrst order of m term y{r)smaii- But the second 
eigennumber called the flrst kind of hidden radial quantum number only appears in the flrst 
order of the m term. As we neglect small mass m, the primary radial quantum number 
only appears in wave function. However, when we include the small mass m, the second 
eigennumber, 1^* kind of hidden radial quantum number is created. When one sees any 
other special function, for example, Laguerre and Associated laguerre function, Legendre 
and associated Legendre function, hypergeometric function, Kummer function, etc, those 
functions only have one eigennumber. This is quiet a strange function. Also, surprisingly, 
the story does not end here! In the new special function, there are two terms which are 0*'' 
order of m term, y{r)dominan. and 1"** order of m term, y{r)smaii- Actually, higher order of 
m terms do exist. But the mass m is extremely small. So We neglect the m terms that are 
more than T^"^ order. We only count it up to a 1^* order of m term. The full description of 
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function y{r) include all higher order of mass m in the following way. 



oo oo 



i=0 n=0 



n=0 



n=0 



n=0 



n=0 



(26) 



If the function y{r) of (126|) is infinite series, then the wave function ?/'(r, 5, 0) will blow up 
as we see in ( 1T4|) . So all of each of summation in (|26|) must be a polynomial. Then, ( 126|) 
become 



No 



N2 



i=0 n=0 n=0 n=0 ra=0 n=0 

(27) 

Then, for example, all possible eigenvalues and eigennumbers of m*'* order term up to forth 
order of m are 



mO order term; j^o =46^/ + 2no - |) tiq = 1, 2, 3, 4, ■ ■ ■ | 



(28a) 



m^ order term; 



El 



El 



46(^/ + 2no- l) no = 1,2,3,4, 
46(^/ + 2ni + i) ni = 1,2,3,4, 



(28b) 



m order term; < 



^2 
-^0 



46(^/ + 2r2o- no = 1,2,3,4, 
46(^/ + 2ni + m = 1,2,3,4, 
46(/ + 2n2 + 



= 1,2,3,4, 



(28c) 



m^ order term; < 



El 



El 



El 



El 



Ah{l + 2no-\) no = 1,2,3,4, 
46(^/ + 2ni + i) m = 1,2,3,4, 
46(^/ + 2n2 + 
46(^/ + 2n3 + 



n2 = 1,2,3,4, 
713 = 1,2,3,4, 



(28d) 
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where 



m order term; < 



Eq 


= 4:0 \ 


I + ZUq — 2 1 


no = 


1 O Q /I 

1,2,6,4:, 


El 


A h 
= 40 


I + Zrii + 2 j 


rii = 


1 O Q /I 


< El 


A 1 

= 4o 


7 1 o 1 3 \ 
/ + 2^2 + f 1 


n2 = 


1,2,3,4, 


El 


= 4bi 


+ 2^3 + l) 


n3 = 


1,2,3,4, 


El 


= 4bi 


^/ + 2n4 + 


rii = 


1,2,3,4, 



nj < Uj where i < j and i,j = 0, 1, 2, • • • 
primary radial eigenvalue 



Ei 



E! 
no 



type hidden radial eigenvalue 



primary radial quantum number 



type hidden radial quantum number 



(28e) 



(29) 



We see z*'* term of m has (i+1) different eigenvalues from f l28al) - fl28el) . Actually, we don't 
have to think about eigenfunctions of all higher order of m^^ term, since these functions 
create extremely small vibrations and variation because of small mass m. So zeroth order of 
m and 1st order of m terms are sufficient. However, we obtain many eigenvalues whenever 
the order of m*'^ term increases. These are reasons why there are many linearly increasing 
lines with same slopes including bunch of eigenvalues corresponding to fermions and bosons. 
From fl28al) - (l29l) . we can obtain quiet exact mass formulas of many elementary particles. 



II. ASYMPTOTIC BEHAVIOR OF GRAND CONFLUENT 
HYPERGEOMETRIC FUNCTION 

Now let's generalize this new special function. Suppose that there is a second order 
differential equation which is 



x'^y (x) + aoxy (x) + {aix'^ + bix^ + c\x^ + d\)y{x) = where , 



do, CLi, hi, Ci, (ii = 3? 
< a; < oo 



(30) 
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(130|) is equivalent to ([T]). All coefficients in the above exactly correspond in the following 
way. 

ao < — > 2 
61 < — > —mb 



Ci i — 7- y-^ m 

di i — > -l{l + 1) 



X — f r 



When a function y{x) = e5v^('^'+^) (_(ao-i)±V(«o-i)2-4di)^^^^j ^^^^^ ^^^^^ 

x(7 (x) + (/ix^ + ex + v)g {x) + (i7a; + eu)g{x) = (32) 

where, 



v=\± V(ao - 1)2 -4rfi (33) 
6^ 

= ci - + «vV(2 ± v/(ao - 1)2 - 4rfi) 
u; = i(l±v/(ao-l)2-4cii) 

( l32l) is equivalent to ([6]). All coefficients in the above exactly correspond to the following: 

yu i — > -b 
6 i — > —2m 
ui — ^2(/ + l) 

n^Q = b(i + - 



4 V 2 

UJ < — > (/ + 1) 



X 
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We first suggest that ||| = using the function g{x) as Frobinous series in 



I2l 12^ 

f l32|) . we have two indicial roots, Ai = and X2 = 1 — jy. And, recurrence formula for all n is 



Kn — An + 



Br, 



: n > 1 < 



An 
I ' = 

Qi. = 

Co 



- Cn 
Cn-1 

£{n+\+Uj) 

(n+A+l)(?i+A+v) 

r2+/i(n+A-l) 
{n+\+l){n+\+u) 

e 

'2 



(35) 



Let's test for the convergence of the function g{x). As n goes to infinity, recurrence formula 
is approximately equal to 

Cn+l ^ -^Cn-l (36) 

n 

Now we can describe g{x) as power series by substituting ( 136|) into Frobinous series. For 
simplicity, we suggest Co = 1 



lim(7(x)~ VC2„x2" + VC 



n 



2^* 



2n+l3; 



= 1 + v^e' 2 
> (1 - ^£x)e-^'^"' 



--/.x^ Erf 



exe 2 

2 



(37) 



In the above, Erf 



1 2 

2 ^^tjC 



is error function. Substitute ( 1371) into ?/(x). 



limy(a;) ~ a;K-("«-i)±V(-o-i)^-4'ii) gtv^x2 _ 



z6i 



a;e 2 



+ (z)=^/2^(ai)i/^xe-5v^" Erf 

^ g-f >rx2^i (-(ao-l)±v/(ao-l)2-4di) 

(I) As ai = real positive 

(a) If |( - (ao - 1) ± v/(«o-l)'-4rfi) < -1, 
As n > 1 and in (IHl 



ibi 



zX 



(38) 



ihi 



hm y(x) > x^ (-K-i)±v/K^Tj^) 1^ 1 



X — 7- 00 



(39) 
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As X of Erf 



(i)3/2(ai)V4a; 



function goes to oo, it's imaginary part oscillates around zero. 
As ai increases, its oscillation are really close to zero, and it's real part oscillates at around 
— 1. Also, as ai increases, its oscillation is really close to —1. As x goes to oo in (138|) 

ibi 



lim y(x) ~ a;l(-(«o-i)±V(-o-iP-4rfi) ( gfv^x^ _ 

n»l ^ \^ 



2^07 



+ (i)3/27^(ai)^/Se-^v^" Erf 







(40) 



(b) If 1( - (ao - 1) ± ^{ao - 1)2 - 4di) = -1, 
turn out to be 



\imy(x) ~ ie-5v^"' - ^^g-^v^"' + (z)3/2y^(ai)^/'e-5v^"'Erf 



(^)3/2(ai)i/^a; 



> e 2 



1 ibi 



As X-)- in (jH]), it then yields 



lim y{x) > e 

n>l 



00 



(41) 



(42) 



As X goes to zero, the function y{x) become divergent. Also, as x— )■ 00 in pij) . suggested 



by III = I ^^7^ I <^ 1, it is simply approximately 



lim y{x) ~ (2)3/2y^(ai)i/^e-iv^"'Erf 



n>l 



fz)3/2(ai)i/4x 



(43) 



(c) If -1< 1( - (ao - 1) ± v/(ao-l)'-4di) < 0, 
As X goes to zero and 00 in 



lim 



y{x) > e-§v^^'a;K-('^«-i)±V("«-i)'-4di) ^ 



ibi 



X — )■ 00 



(44) 



(d) If K - (ao - 1) ± v/(ao - 1)=^ - 4rfi) = 
(l38|) turns out to be 

ibi 



limv(x) ~ 62^^""^ — — 



2;e-5v^-' + (i)3/20F(ai)V4a.e-5v^-'Erf 



Uf^(aiy/^x 



- i . /7TT,.2 

> e 2 



^^^■^ f 1 _ 

2^ 



(45) 



As X goes zero and 00 in (145 



lim ?/(x) ~ 1 

n>l 



lim y{x) —7- 00 

x— >-oo 



(46) 
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(e) If i( - (ao - 1) ± v/(ao-l)'-4rfi) > 0, 
As X goes to oo and zero in (!38l) . 



lim y{x) > e 



X 



l(_(ao-l)±^(ao-l)2-4di) A _ J^^^A 

V 2^ / 



— 7- OO 



(47) 



z6i 



+ («)3/2yi(ai)^/^xe-^v^"'Erf [(z)3/2(ai)i/4x' 



^ 



(48) 



(II) As ai = 

(|38|) simply turns to be 



lim y(x) ~ a;5 (-('^o-i)±VH-i)^-4rfi) i i _ 



(49) 



We suggest that 



^ 1. Then as ai = in the second term of the bracket in ( H9l) . 



—7- OO. The function y{x) will be divergent no matter what the value of x is. 



I 2 I I 2y^l 

Therefore, there are no any independent solutions at all in the case of ai = 0. 
(Ill) As Oi = real negative 
Plug ai = — |ai| into fl38|) . 

\imy{x) ~ a;l(-("o-i)±VK-i)2-4di) ^-^^\x^ ^ 0F(|ai|)i/^a;e5^"' Erf k|ai|)^/^x' 

n>l Y L 



2v^ai| 

> glA/KI^'^,! (-(ao-l)±V(«0-l)2-4di) (^-^ 



&1 



(a) If K - (ao - 1) ± v/(«o - 1)2 - 4rfi) < 0, 

As X goes to and oo in ( 150|) . the function y{x) is divergent. 

»i ^ ^ V 2v^ ^ 



(50) 



lim 

n>l 
x—^oo 



— i- OO 



(51) 



(b) If i( - (ao - 1) ± ^{ao-iy-M^) = 0, 
As X goes to and oo in ( l50l) . 



lim y{x) ^ 1 



lim ?/(x) —7- OO 

ra>l 



(52) 
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(c) If i( - (ao - 1) ± v/(ao-l)'-4rfi) > 0, 
As X goes to and oo in (150|1 . 



\imy{x) — )■ lim y{x) — )■ oo (53) 

n>l n>l 

We check all possible tests to determine if an infinite series of function y{x) converges 
or diverges on the above. Now let's test for convergence as the function y{x) as g{x) is 

N 

polynomial fot the case of ai =real negative. Substitute g{x) = ^^C„a;"^'^ into y{x). 

n=0 

N 



(-{ao-l)±^(ao-l)2-4di) ^ (^7^3;"+^ (54) 



y\x) n 

n=0 



(a) If 1( - (ao - 1) ± v/(«o-l)'-4c/i) < 0, 
As X goes to and oo in fl54|) . 



lim ?/(a;) — )■ oo lim ?/(a;) — > (55) 

n>l ?i>l 



x-s>0 



(b) If i( - (ao - 1) ± v/(«o-l)'-4rfi) = 0, 
As X goes to and oo in ([5i 



lim?/(x) — )■ 1 lim y{x) — > (56) 

n>l n>l 
X— >0 X— >oo 



c) If i( - (ao - 1) ± v/(ao - 1)' - 4di) > 0, 
As X goes to and oo in fl51|) . 

lim y{x) (57) 

n>l 
x-s>0 

X— >oo 

Now let's try to investigate independent solutions of a function g{x) for two cases which are 
polynomial and infinite series. First of all, if we choose boundary conditions that a function 
g{x) goes to zero as x goes to oo, and it is convergent as x goes to zero. Then, the function 
g{x) must be polynomial. The boundary condition for the polynomial is the following: 

limg{x) — convergent 

(58) 

lim g{x) — 7- 

X— >-oo 
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and the necessary conditions of it is 



fli = real negative 

|( - (ao - 1) ± y/{ao-iy-Adi) > 
< 1 



(59) 



I I 



After we develop independent solutions for the polynomial, we can expand them as infinite 
series in simple ways. But even if this process is simple, it does not mean that polynomial 
series is just the special case of the infinite series. They are independent each other. When 
you try to find the solution of any differential equation, first of all, you must consider what 
physical circumstance and mathematical condition make solutions as polynomial series and 
or infinite series. As we know, there are two indicial roots which are Ai = and A2 = 1 — z^. 



III. FINITE SERIES 

A. u = non-integer 

1. As Xi = 

a. Power series 
From (I35D 

[n + l){n + u) 
{n + l){n + u) 

Put n=0 in flGOaj) . We get Aq = -^e = -| = gi. From (I35]), we have /Tq = § which is 
equal to Aq. Plug n = 1 into recurrence formula. 

As we see in (|6T|) . An includes the first order of | in which is has an extremely small value. 
Then, we argue that \Ai\ ^ |^| = (I6T]) is approximatively equal to i^i ~ ^ = —j^- 

By using this process, we can simplify f p5|) by giving Kn in terms of An and -B„ instead of 
Kn-i up to the first order of | . By using = -tt^, even and odd terms of C„ are same 
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as (fT6|) . We can describe g{x) as power series by using ( !T6|) . 

oo oo 

g{x) = g{x)domin. + g{x)smaU = ^ C2nX^'' + ^ C2n+ia;^"^^ 

n=0 n=0 
J- oo n—1 s f 

= Co|l + ^ n ^2^+1^2" [+ CoMoX + {AoB, + As^Ox^ 

n=l fc=0 ^ 
oo p n—1 n—1 y n—1 J — 1 \ n—1 

+ XI ^0 n ^2p+2 + J2 (^2j n ^2p+2 n ^2fe+l ) + ^2n n ^2p+l 

r,—0 L n—n o — 1 V r,— o I fl / r,— O 



(62) 



p=0 j=l ^ p=j fc=0 ^ p=0 

By using similar process from the previous case, there are two eigennumbers which are 
— — = 77-0 — 1 = |ao| where rio = 1, 2, 3, ■ ■ ■ (63a) 

- + ^ j = ni - 1 = |q;i| where ni = 1, 2, 3, ■ • • (63b) 



f l63ap makes i?2fc+i term go to zero at certain value of k where k = 0, 1, 2, ■ ■ ■ . And f l63bp 



makes -82^+2 term go to zero at certain value of k. fl63ap and (]63bp make g{x) function 
a polynomial series. The g{x) small term is extremely small value relatively compared to 
9 (x) domin. because it includes An term having |. First of all, add fl60ap and f l60bp into the 



first term of g{x) in (!62|) . putting Cq = ^^^^rr^' ~ + ^) ^ ~ —hfix'^. 



f ^ 17 f \ r(|ao| +7) (-|"o|)n n (aA\ 

g{x)domin. = i^aol (7; z) = ^^1^7)7^ ^ ^ 



El]) is same as the first kind of confluent hypergeometric function. Substitute f l60ap and 

r(7) ' ' 2> 



(I60bp into the second term of g{x) in fl62|) . plugging Cq = 7 = ^(1 + z/) and 
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. . __r(7±M J V — (-ir(i)K7 - i)!(l«i|)! 



r(7) Xt'o 2(|)(7 - i) + |)!(7 + n- i)!(|ai| - n)! 

<n + ^) (-l)"(7-l)!(|ao|)! 
^2(n + i)(n + 7 - i) (n)!(7 + n- l)!(|ao| - ^)! 

e{j + I) (-1)"(7 - l)!(|c^o|)!(j + |)!(J + 7 - |)!(ki| - j)! ' 

t^^2(j + |)(j + 7 - i) (J)!(7 + J - l)!(l«o| - j)!(^ + + 7 - i)!(l«i| - n)\^ 



(65) 



17 



Putting |ao| = 0, 1, 2, ■ ■ ■ in (165|) . we yield 

^ (n + |)(n-i)!(|ao|)!(|«o|+7-i)! 

fi- (2;) sman = -X^}^ 



„=o + ^ - i)H!(-|)!(l«o| - n)\{n + 7 - 1)! 

( -l)^(|)!(n + 7-|)!(|ai|-n)! 
(k + ma^\ - k)\ik + ^ - l)\ 



(A: + i)!(|a,|-fc)!(A: + 7-i)' ^''^ 
As we see in f l66|) . maximum value of index n is |q;o|. The range of index k is < A; < |q;i|. 
In other words, < n < \ao\ < k < \ai\. Then we obtain |ao| < If |cto| > the 

function g{x) will be infinite series. Then, the function y{x) will blow up as we plug g{x) 
into it no matter what the value of x is. Such solution does not exist. When we see the 
second summation of flBBl) . we can shift index k to zero at the beginning of summation. Then 
we can replace the interval of index k by < A; < |q;i| — fl66l) which is simply described as 

_ ^^ r(|c.o|+7) ^ tMk.n ' v^" in + |)r(n + i)r(n + 7 - - 
2 r(7) ^ -!(7)n ' r(A; + n+|)r(A; + n + 7 + i) 

(67) 

|ao| 

We see the first summation of f l67p which is ^(I°o|+'t) z"-, it is the first kind of con- 

n=0 

fluent hypergeometric polynomial of degree |ao| which is denoted as F^aoiij] z) surprisingly. 



Substitute (EID and (ETD into ([62]). 

g{x) = Q>V|„o|,|„i|(^|ao| = "-0 - 1, l^il = ^1 - 1,7 = ^(1 + = ~^/^^ 

= ^|ao|(7;^) - 2^11|„^l W'^) only if |ao| < |ai| (68) 

where, 

yrM. . _ r(|«o|+7) ^ (-|ao|)n „. ' v^" + !)r(n + |)r(n + i-l)in- \a,\), , 
AAKr^'"'^ r(7) A. ^!(^)^ ^ r(A; + n+|)r(fc + n + 7 + i) 

(69b) 

(jHE]) denoted as QW|q,q|_|q,^| ^|ao| = no — l,\ai\ = rii — 1, •y = ^{1 + u); z = — is called 
as the first kind of Grand confluent hypergeometric polynomial of degree |ao| and |ai| with 
the first order of |. 

We also we obtain two eigenvalues from f l63ap and (163bp . which are 

Ql = -2fi\ao\ where |ao| =0,1,2, ■■■ (70a) 

nl = -2fi(^\ai\ + ^^ where |ai I = 0,1,2, ■■■ (70b) 
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b. Integral formalism 

The solution of the Laguerre differential equation is 

A:=0 

And the solution of the associated Laguerre differential equation is 

L'(z) V i-iyin + k)\ . ^ e^z~^ ^ 

( I69bl) might be described in the following way: 

I I I I 

n: (^;^)=Ee»'(7;.)nr"(o..) 

n=0 

where, 

r(7) n!(7)„ 

" ^ r(A; + n + |)r(A; + n + 7 + i) 

Plug |ao| = into ([68]), fl69al) and fl73al)- fl73cll 



z' 



z^ 



k=0 

The beta function is 



^^'^^ r(p + g) io Jo (l + tn" ^^"^^ 

And, 

By using (176!) into the second term of (17^ on r.h.s 



A;=0 

Replace p and q by 7 — | and A; + 1 in (!75|) . Substitute it into ( I77|) . 



fc=0 
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Replace n and z by |ai| and A = z{l — t){l — p^) in (17T|) . Then, plug it into (178|) 

nr'(^'^)=^ f f dtt^-H\^,\{\) 

u / J_i Jo A=z{l-t)(l-p2) 

Plug ([79]) into dUD 

/I 
dp / dtr-iL|,,|(A) 
-1 ''0 A=2:(l — t)(l— p2) 

Plug |ao| = 1 into ([MD, and flTgiD- dTScj) 

QWi,|.i| (7; z) = F,{r, z) - |x|fo1(7, z) nl,"' (7; z) + F/(7, JJ^' (7; ^) } 
Replace k by A; + 1 in fITB]) , and plug it into the second term inside of bracket in f lHT]) 



(79) 



(80) 



^1) 



|ai|-l 
A;=0 



;-i)'=r(|ai|)r(7 + i) 



2mA: 



r(|ai| -A;)r(A; + 7 + |)r(A; + 2) 



[z{l-p')] 



(82) 



Replace p and q by 7 + ^ and /c + 1 in (I75|l . Substitute it into fj82|) . 



-lfn\a,\)[z{l-t){l-p^)] 
r(|ai| - A;)r(A; + l)r(A; + 2) 



(83) 



Replace n, k,j and z by |q;i| — 1, 1, k and A = z{l — t){l —p^) in fl72l) . And plug it into (|83|) . 



Substitute (1791) and (IHl into (EI 



QWi,|„,|(7;2:) =Fi(7;2;)-|x|fo1(7,2)| /^^^P^ rft t^-tL|,,|(A) 

(7,^)^(1 + I) j\{l-p') j\tt^-Hl^^_,{\)^ (85) 

As |ao| = 2 in f l68|) . f l69ap and f l73ap - fl73cp . the result of its solution is similar to the previous 
case as |q;o| = and 1. 



2 |. 

QW2,|„,| (7; z) = F^ir, ^) - 1^ E ^nil, zKn + 

n=0 



ujJ{n+l)T{\ai\ + l-n) 



r(i)r(|ai| + l) 



I' dt t^-^-^'^ £^ dp (1 - p^)"L|^,|_„(A) I 



(86) 



A=2(l-t)(l-p2) 
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According to ([HOD, ([85]) and ([86]), fl69i]) and fl69bl) are 



And 



n::;:(7;^)=Ek°»'(7.^)(«+?) 



l"ol 



-|"ol 



;l"ol+l(l 



r(n+|)r(|ai| + l-n) 

r(i)r(|«i| + i) 



n=0 - - V2 



A=2(l-t)(l-p2) 



We know 



r{n + i)r(|ai| + 1 - n) B{n + \, |ai| + 1 - 



r(|)r(|«i| + i) 

Substitute (!73bD and ([89]) into ([88]). 



n:::(-'^-)HE 



-|«o| 



|«o| 



(-ir(|ao|)!(|ao| +7-1)! B{n + 1, + 1 - n) 



- l)!(|ao| - n)\{n + 7 - 1)! |«i| + 1) 

u (-If (|ao|)!(|ao| + 7 - 1)! + |, + 1 - n) 
2 n!(|ao|-ri)!(n + 7-l)! + 1) 



A=2(l-t)(l-p2) 

Replace p and q by n + | and |q;i| + 1 — into fl75l) . Then, plug it into f l9Q]) . 

1 



n::::(-) 



2' 
|ao| 

?i=0 
|ao| 



-ir(|ao|)!(|«o|+7-l)! 



1 

-1 



(n- l)!(|ao| - n)!(n + 7- 1)! 



a; ^-l)»(|ao|)!(|ao|+7-l)! 
+ n!(|ao|-r2)!(n + 7-l)! ^'^'^^ ' ^ >^ Vl-n(^) 



n=0 



A=z(l-t)(l- 



Integral form of Associated Laguerre polynomial is 



1 
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Replace n, m and z by |q;i| — n, n and A in (192|) . Then plug it into ( 19T|) . 



n::::(^-) 



2mBC^,\a,\ + l)Jo Jo 



Alt 

e I , 9x M 



u 



zst{l — p 



1-u 



\ao\ 



n=0 



(-ir(|ao|)!(|ao|+7-l)! 
(n- l)!(|ao| -n)!(n + 7-1)! 



zts{l — 



u 



l"o| 



+ 9E 



n=0 



-ir(|«o|)!(|«o|+7-l)! 
n!(|ao| -n)!(n + 7- 1)! L 



zts(\ — 



{l-u) 



u 



n-1 



{l-u)\ 



(93) 



A=2(l-t)(l-p2) 



Replace z and A by = — p )jYZuj ^^"^ '^(-'- ~ ~ P ) i^^^^ (03 



n!:::(--) 



1 



27ri5(i, |ai| + 1) 



X (b du 



/ cis s-^(l + s)-(l"^l+i) / rftt^-i / 
io io J-1 



(94) 



u 



can be described as various integral forms of several different special function in the 
following way: 

n'"'(7;^) = ^i^;^^^^ r rrf..-ni + s)-(M+i) [\tt^-l Cdpi , J" , 



(95) 



n!::!(-) 



(2^)25(1, + 1)70 ^ ^ Jo 



X 



dv 



du 



zstvu 



where, 1^3 



y|ao|+i(i _ yyi J mI"iI+i(1 - n) [ (1 - v){l - ti) 

x(iM(l,^,^3)-M(l,|^3))+|M(l,|^.3) 

stv 



zu 



{l-u) 



(1-t) 



(96) 



n|::|(-) 



(l«o|)! 



(27ri)25(|, |ai| + l) 



/ s-^(i + s)-(i"ii+tM dp (i-p'^y-^ (b 

Jo Jo J 



X 



zup 



-zspvu 



due 

mI"iI+i(1 -m) I (7 + i)(l - v){l - u) 
3TT^(7 - + ^,^4)} where, W4 



^(7+ 2' 7 +2' ^4) 



dv 



2(7 -i)' 



(97) 
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X { - 2(i-.)(i-,) ^(2'-'"^' + ^'"^) + 2^(2'-|«il'-^)| 

where, = — -(1 - p ) (98) 

{l-v){l-u) 

In the above, M(a, 6, z) is the first kind of Kummer function which is 

(a) 



M(a,b,z) = V -y^^" = e"M(6 - a, 6, -2) 



n=0 



where, Re(6) > Re(a) > (99) 
And f/(a, 6, z) is the second kind of Kummer function which is 

^' ^) = ^' ^) + ^?^^'"'Af (a - 6 + 1, 2 - 6, z) 

1 (a — + 1) 1 (a) 

= z^'^U{l + a-h,2-h,z) 

_ 1 

~ f(a 

Substitute integral form of the first kind of confiuent hypergeometric function into ([9 



-ztja— l/i I a— 1 



rft e-'H^-'il + 1)"-"-^ where, Re(a) > (100) 







2' 

1 /•! 

3 



X / dt t^-S / dp (t I 1,,^^ (101) 

^ / _ ZStUv{l-p^) C^l 



_ (l-u) 

d°il+i(l - n) \ (1-m)(1-i;) 2 
Description we gave in (p5|) -( l98|) as integral formalism is exactly equivalent (llOip . And 

we obtain the integral representation of the first kind of grand confiuent hypergeometric 

function according to ( IHTll and fllOll) . 

c. Generating function 

Now, let's try to get the generating function of the 1^^ kind of grand confluent hypergeo- 
metric function. First of all, we multiply + 1, -) vj^^^ on both sides of integral 

|ai| = |aol 

form of 1^^ kind of grand confluent hypergeometric function in Table 2 where \vo\ < 1. 



°° 1 11 

J2 B{\a,\ + 1, -) t;i°^'QW|„„|,|„,| (7 = -(1 + uy,z = --fix') = F\^,\{z)l - 

(102) 



|ai|=|ao| 
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where, 



|cn|=l"ol 



11= 5: B(b.i+i.5)"i°''nli;l (■'--) (^o^^) 

\ai\ = \ao\ 

The 1^^ kind of hypergeometric function is 



A;=0 



only if 1^1 < 1, Re(c) > Re(6) > (104) 
Substitute (fTOD into (fT03al) . 

i {\aQ\ + 2 j Jo 

(105) 

Substitute ffTOTD into ffT03bl) . 



dv 



t;l"ol+i(l - 
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Ml + --^oV (l + s-^;o)(l-^;) + 2j/" (^°^^ 
Multiply V ^ with both sides of ffT02D where < 1. 

|ao|=0|ai| = |ao| ^' 

= A - |xB (107) 

where, 



oo |qo| 



|ao|=0 

oo Jaol 
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Plug ffT05|) into ffTOHajl . 



|ao|=0 







= f i^ + n + rn- m-zvoV.nvoVrr f\t ^^-{1 - tyhl - v.t^ (109) 
n!m!(7 + n-l)! Jq 

n,m=0 ^ ' " 

Replace a,b,c and zbyl, n + m + 1, n + m + | and Vq into (11041) . And substitute it into 

(iinH!). 

^ i? (n + m + 1, 1) (7 + n)^{-zv^v^r[v^v^T ^ ^ ^3 

A= > — ^ — 1 — I 2-ri l,7^ + m + l,n + m+-,t;o 

^-^ nl ml V 2 



n,m=0 



oo 



V- ^(n + m + J + 1,1) (7 + n)ra{-z) ^^^^^^j 

n,m,j=0 

We can describe f lllOp in a different way. First of all, there is the 1st kind of Appell 
hypergeometric function which is 

/3, 7; x, y) = i^Mp^^^./ for < 1, \y\ < 1 (111) 

The function Fi can be expressed by simple integral 

r(7) Jo 

where Re(a) > 0, Re(7 - a) > (112) 



Put a = n + l, 7 = n + |, /3 = 1, /3' = 7 + ri, x = t'o and y = VqVi in (I112p . And substitute 
it into ffM . 

. ^r(i)(-;2t;ot;i)"' , 3 , 

— r( _L 3^ 1 (ra + 1; 1, 7 + ra; n + -; wq, ^^o^^i j 

n=0 ^ 
_ ^(^ + + J + 1. I) (7 + .^^o^ 

n,j,fe=0 



(11131) is exactly equivalent to the second line of (IllOp . Substitute (I106p into (llOSbp 



B = T,{s,p,t) + T2{s,p,t) (114) 
25 



where 



ri{s,p,t)= ^^^'^"p^i ^ dsskl + s)-^l + s-vo)-' dtt^- 



2 



X / dpp{l-p)—2e '-^^} (115) 



cj poo pi 

T2{s,p,t) = —^— dss-^l + s)-kl + s~Vo)-' dtt^- 

(1 - Vip Jo Jo 



X / dp{l-p)--2e ^-"iJ (116) 



By using the first kind of Kummer function, 



j\pp{l-prh~"P = e-"|2M(i I a 



2,^/3 5 
3 V2'2' 



ZVq f tsv 

wnere a = — ' ' 

14 

Plug (dUD into (fnSD . 



where a = + ^| (117) 



T^{s,p,t) = -fl{s,t)+-f2{s,t) (118) 
^^(s,t) = jf^^^^^p— [ ds s^l + s)-^l + s -Vo)~'^e~^^ (119) 



(1 -fi)^+i 



tsvi 

2' 2' 1 + s-vo l'' ' 1 



72(^, t) = 3(i^r°[')Vi ds si(l + .)-! (1 + . - t;o)- (120) 

X r t-^e^ (^-^)'Mf ^, ^, {(1 - t) + 

There are addition theorem for M{a, b,x + y). 

M(a, b,x + y) = y^ M(a + n,b + n,x) (121a) 

^ (o)„n! 



n=0 



M(a,6,x + i/) = e^^^^-^|^^^-^M(a,6 + n,x) (121b) 

Put a = |, 6 = f, X = and y = -i^(l - into ([mbD, and plug it into 
(11201) and replace s by j — 1 in it. 

^ (7 + ^ + i)(i)n io V y 
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Put a = I, 6 = f, X = j:^^ and y = -t^(i - T^)t into ([mbD, and plug it into 
(11201) and replace s by j — 1 in it. 

^ 2 g (,_,)-.^...(^ _ 1 

3^ (7 + ^ + l)(|)n io ^ ^1 

5 ztvo 

Vot 

We know that 



m(„ + 1,„ + I ^) = V (" + _ .„„()-". (124a) 



2 1-votJ ^ (n + |)mm! 



m(„ + l.n + 5. = f (!i±iUz£!^r(l - „„r" (124b) 



Plug ffT2lil) and ffT24bD into([l22]) and ffT23|) . And Substitute ffTT2|) . putting a = m + 1, 
/3 = n + m + 2, /3' = — n, 7 = m + |, x = fo and y = ^ into them . 

^ " ' nho (7 + - + i)(|)> + iUr(m + |) 

/ 5 1 \ 

xFifm + l;n + m + 2, -n;m+ -;fo,— j (125) 

2 ^ (-l)"+"^+^(l-t;i)-("+^+^)(-2^o)"+'"+'(wi)"+'r(|) 

^'^''^^ ^ ~3^^^ (7 + n+i)(|).(n + |)„r(m+|) 

/ 5 1 \ 

xFi ( m + 1; n + m + 2, — n; m + -; fo, — ) (126) 
V 2 vi/ 

Substitute dHSD and dm} into (fTTS]) . 

A_ V" (-i)"^'"'''(^ + ^ + i)r(|)r(|)(i-t;i)-("+^+'n^^o)"+'"+'K)"+' 

^ ~ nho (7 + n + i)r(m + |)r(n + m + |) 

/ 5 1 \ 

xFWm + 1;?2 + m + 2, -n;m + -;t;o, — j (127) 

Replace s and p by ^ — 1 and 1 — in (11161) . 

12(5, p, t) = -^^^-y y dv {l-v)-2{l-VoV)-^ J dtt^-h '-'^""V ''(^-"i) 



(1 - io 



dv {l-v)-^{l-vov)-^ dtt^-^ 
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Put a = 1, 6 = |, X = and y = j^jl - i^pz^ji^ into (fmil) . and plug it into 



^ . ^ ^ vl{l - vi)-^{-zv^Y+^{n + 1) 



(1 - ^i)^ (7 + ^ - |)(|)n(n + \)m m\ 

dvv'''{l-v)-'^{l-Vov)-^''+"'+^'>(^l- (129) 



Put Q; = m + 1, 7 = m + |, /3' = — n, /3 = n + m + l, x = fo and y = in f lll2p . Then plug 
it into ffT29D . 



>^ (-l)"+"^r(|)(n + 1)^(1 - t;i)-("+^)(^t;o)"+"^(t;i)" 
- c^^Z.^ ^ ^ _ i)(3)^(^ + 3)^r(m + f ) 

xFi(^m + l-,n + m + l,-n-,m + ^-,Vo,^^ (130) 



Substitute (fT27D and (fT30|) into ffTTij) . 



oo^ (_l)n+m+i(^ + ^ + l)r(i)r(|)(l-i;i)-("+7+i)(z^;o)-+-+i(i;i)"+i 

^~„^o (7 + ri + i)r(m + |)r(n + m + |) 

/ 5 1 

xFi (m + 1; n + m + 2, —n; m H — ; Vq, — 
V 2 v^ 

oo 



(-l)"+™r(i)(n + 1)^(1 - vi)-(«+^)(zwo)"+'"(^^i)' 



„,„=o (7 + ^ - DilUn + l)mT{m + |) 

xFi(m + l]n + m+l,-n-,m + ^-,vo,^^ (131) 



Plug (I113p and ( I13ip into (I107p . Then we obtain the generating function for 1^^ kind of 
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grand confluent hypergeometric function: 



(|ao|)! 

— r( M 3\ ^1 (n + 1; 1, 7 + n; n + -; Wo, VqVi ) 

n=0 U^+2i ^ 



-X 



g + + l)r(i)r(|)(l - i;i)-("+T+i)(zt;o)"+'"+iK)"+' 



2 [^^^ (^ + ^+ i)r(m + |)r(n + m + |) 

/ 5 1 N 

xFi m + 1; n + m + 2, — n; m H — ; Vq, — 
V 2 Vi^ 

(_l)n+-r(l)(n + 1)„(1 - t;i)-("+^)(^t;o)"+'"(^i)' 



„,^=o (7 + n - i)(|)„(n + |)„r(m + I) 

/ 3 1 \ 

xFWm+l; n + m + 1, -n.; m + -; fo, — 
V 2 fi/ 

_ + fc + J + 1, |)(7 + r(i)r(|) - 

nj,A;=0 n,m,j,k,s=0 
^_l^n+m+l^^ ^ ^ ^ ^ ^ ^ ^)^(^ ^ ^ ^ ^ 2),(-n)fc ^"+'"+^ ^n+n.+, + l^,n+.-fc+l 

s! j! fc! (7 + n + |)(m + |)!(n + m + |)!(m + |)j+fc 
, ^ (-l)"+-(n + 7)s(r^ + l)^(m + l),+,{n + m + l),{-n), 

.!j!A;!r(|)(7 + n-i)(m + l)!(|)„(n + |V(m+|U. / ^'^^^ 

We can describe tlie generating function of it as the integral formalism by using fllOQp . flllSp 
and f lTTHl) . 



f f (ill)'' "^'"^'Q>^l°.M°■l (7 - -l^- 

=0 |oi| = |cio| 

dt (1 - - voty\l - vovitye'^^^' 

du [1 — u)~^{l — Vqu)'^ / dtf~^^ / p)"^ 



e X 



^0 JO 

""0 J ^ f(n-ui) 



'2 (i-t,i)(i-t,t,i)/ ^ ' 



d. Orthogonal relation 

From the differential equations satisfled by QW|Q,o|ja^|(-2) and QV^ \i3oi\i3^\{,z) , namely, 

^2>V|„o|,|„,|(^) + + 5X + z/)QW;„„|,|„^|(z) + + £u;)QW|„„|,|„,|(2;) = (134a) 
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n 



\ao\,\ai\ 



^Q^lMmi^) + (/^^' + ^^ + ^)Q^'\Pol\M^^) + i^\M,\M^ + ^^)QW|/3oM/3i|(^) = (134b) 
where, 

multiplying fll34ap and (1134bp by QyV\f3o\,\i3i\{^) Q^\ao\,\ai\{z) respectively and subtract- 
ing, we have 

dM{z) 



X- 



dx 



+ (/ix +ex + u)M{z) = mpo\,\M - ^l«o|,|ai|)a;QW|«o|,|«,|(z)QW|^o|,|/3i|(^) (136) 



where 



M{z) = QW;„„|,|„^|(^)QW|^„|,|;,,|(z) - QW;^„|,|^^|(^)QW|„„|,|,,|(; 



(137) 



Multiply x" ^62^^^ on both side of (11361) Then integrate it with respect to x from to 
oo. 



(^l/3ol,l/3i| - ^l«ol,|ai|) / ^'"Q>^|ao|,|«i|(^)2>V|/3oM/3i|(^) 

Jo 



>oi,|ai|(^)2^l/3oi,i/3i|(^) - 2>^I/3o|,|/3i|(^)2>Vn|,K|(^)) 
Therefore, only if |q;o| 7^ |/3o| and |ai| 7^ then 

dx x^e^'^^-'+^^QW|„,|,|„,|(z)QW|^„l,l^,l(z) = 



(138) 



(139) 



Let's think about the case in which are |ao| = \Po\ and = First of all, from (!68l) 



1 2 



^Nl(^) 



(140) 



e is extremely small. So anything beyond the second order of e is negligible in f ll40p . Multiply 



j^u^lfix^+ex f|X4Qp and integrate them with respect to a; = [0, 00] 



dx a;^e5^^'+=^ 



1 2 



^1 - 



where 



El 



dx a;^e^'^^'+^^ 



^|«o|(^) 

dxx'^^'e'^'^^'^^^F^^.^iz) n|^;| (^) 
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(141) 

(142a) 
(142b) 



Apply both side of generating function of the first kind of confluent hypergeometric function 
into ffT42il) and ffT42bl) . 



27-1 (_l)Nl27-ir(|ao| + 7- i)r(|ao| + 7 + 1^ 
E,^^T{\ao\ + m\ao\+7)+s ^ ^ (143) 



r^r(7-i) 



7+ 



|ao| 



2-'-t r(|a„| + 7) (-|aol)„ (" + |)r(n + i)r(n + 7 - - |ai|)t 



|ai|-n 

r(7) ^ -1(7)7 h ■ 



r(A; + n + |)r(A; + n + 7 + i; 



;_l)l"olr(n + A; + |)r(n + A; + 7 + i) 
r(n + /c- |ao| + 1) 
2 (-l)l"«lr(n + A; + 7 + l)r(n + A; + 2) 



(144) 



^/2H r(n + A;- |«o| + 2) 

We neglect any terms which include more than second order of e, extremely small, in fll43p 
and fll44p . Substitute fll43p and f ll44p into fll4ip with neglecting any terms larger than 
second order of e. Therefore,the orthogonal relation of it is 



97-1 f_lV°ol27-i fr( |ao| +7 - i )r( |ao| +7 + I ) 

'-r(KI + i)r(KI+7) + ^-^ i) ^ J V V 



r 7 



(145) 



l"o| 



r(K|+7)^(-|Qol) 



E 



|oi|-n 



n + f )r(n + i)r(n + 7 - i)(n - |ai| 



..0 2T{k + n + l-\ao\) 



r(7) ^ -!(7)n 

Also, if there is an analytic function \l/(x) having first and second continuous derivatives in 
[0, 00] and approaching zero when a; — )• 00, they can be expanded in terms of QyV\ao\,\ai\{l', z) 
by using fll45p : 



CO 00 11 



fix 



\ao\=0 |ai| = |aol 



where. 



A 



|ao|,|ai| 



dx x'^e5^^'+-Q>V|,„|,|,^|(7;z)vI/(x) 

[ r(|ao| + 7-|)r(|«o| + 7 + i) 
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-1 



-r(|«ol + i)r(|«ol+7) + ^- 



(146) 



.l)hol27-^ 



r(7-|) 



r(|«ol+7) 



hoi 

E 



•|ao|)r 



r(7) t'o '^'(^)- fc=o 



|ai|-n 

E 



n + f )r(n + l)T{n + 7 - |)(n - 
2r(A; + n + I - |ao|) 



-1 



(147) 
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2. As \2 = l-v 



a. Power series and integral formalism 

From the previous case, we get the first kind of independent solution of grand confluent 
hypergeometric function, putting Ai = 0. As we put \2 = 1 — v into 

A=i-. = \ 148a 

[n + \)[n + 2 — v) 

[n + l)[n + 1 — v) 
And there are eigennumbers where 7 = |(1 + z/) which are 

-^ + 7-1 = ^0 whereto = 1,2,3,- ■• (149a) 

-^ + 7-^ = ^1 where^i = 1,2,3,--- (149b) 

From the above two eigennumbers, we obtain two eigenvalues. 

Ql = -2fi{^o + 1-7) where = 0, 1, 2, ■ ■ • (150a) 

nl = -2fi (^1 + ^ - 7) where t/-! = 0, 1, 2, ■ ■ ■ (150b) 

Substitute fll48ap and (]148bp into (!62|) . Then we obtain the second independent solution of 
grand confluent hypergeometric function by using similar process as we did before. 

7^W^o,^l (7 = ^(1 + ^); ^ = -^/ix^) = ^'-^|a^„(7; z) - /\]^; (7; z) I (151) 



where. 



r(^o + 2 - 7) >fi (-v^o)n ^„ 
r(2-7) h^i^^ ^ ^ 

* . , _ r(V;o + 2-7) (-^o)„ ^„ 

/Vo^^' ^ r(2-7) ^^!(2-7)n 

^ t^" (I + n + 1 - 7)r(n + |)r(n + | - 7)(n - V^i). ^^, 

^ r(fc + n+ |)r(A; + n+ I -7) ^ ' 



k=0 

And its integral representation is 



(154) 
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where, 



a::(-) 



(27rz)2fi(i ^1 + 1) 



noo pi pi 

Jo Jo J- 



dp 



X 



dv 



du 



2 ■ 



zstuv{l — p^) 
{l-u){l-v) 



(155) 



f ll55p also can be described as various integral forms with several different special functions 
as we did before in the case of Ai = 0. Also we can obtain generating function and orthogonal 
relation of 2^^*"^ kind of grand confluent hypergeometric function by using similar process as 
we got them of 1^^ kind of grand confluent hypergeometric function: look at (11321) . (11331) . 
(11461) and (I147p . We skip these parts because of huge details. 



B. 



integer 



When V is integer, one of two solution of the grand confluent hypergeometric equation 
does not have any meaning, because y2{x) = T^Wvo.V'i ~ + h');z = — can be 
described as QW|oo|jai| ^7 = |(1 + u); z = — i/ix^j as long as |Ai — A2I = — 1| = integer. 
If z/ = 1, the two solutions are identical. Then two solutions can be solved by Case B) of 
theorem 1 on the below, and if u is any integer except 1, then it can be calculated by Case 
C) and Case D) of theorem 1 in the appendix. Now let's investigate the second independent 
solution as u being an integer. 

(i) As u = negative integer, 0, and 1 



Theorem 1 P(x)0 + Q{x)^ + R{x)G = 

Suppose is a regular singular point ; G{x, A) = Yl'^=o Cn{^)x^^^ 
(l)As Ai < A2.' distinct roots of the indicial equation. 

Case A) If Xi — X2 is not an integer, there are 2 independent Frobenius solutions 



gi{x) = G{x,X) 
Case B) If Xi - X2 

Case C) IfXi 
9i{x) 



A=Ai 



G{x,X) 
X2 is ar, 
G{x,X) 



g2{x) = G{x,X) 
0, there is a Frobenius solution 
92{x) - 



\=\2 



dX 



A=Ai 

A2 is an integer, 



A=Ai 



A=Ai 



92[X) 



dX 



A=A2 
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(2) As A2 < Ai; distinct roots of the indicial equation. 

Case D) If X2 — Xi is not an integer, 



giix) = G{x,\) 



A=Ai 



g^(x) = G{X,X) 

A=A2 

Case E) If X2 — Ai is an integer, 

/ X _ 0((A-Ai)G(:r,A)) 

92[x) 

A=A2 



5-1(0;) = ^(a;, A) 



A=Ai 

Since u is smaller than 1, A2 = 1 — is positive integer. So Ai = is smaller than 
A2. By using the Case C) in Theorem 1, gi{x) = G(x)|a=Ai=o is equal to QW|q,q|,|q,i| ^7 = 
|(1 + u); z = — Also, since z/ is equal to 1, Ai is equal to A2. So by using the Case 
B) in Theorem 1, gi{x) = G(x)|a=Ai=o is equal to Q>V|a(,|,|ai| ^7 = |(l + z/);z = — i/ix^j. It 
is same as before. Only difference is that in the previous case u was equal to or negative 
integer, and in this case u is equal to 1. Now, let assume that the second independent 
solution is 

g2{x) = u{x)Q'W\ao\,\a^\ (l = ^(1 + i^)] z = ^^f^x^) (156) 
Substitute f ll56p into the grand Confluent hypergeometric differential equation 

X {u{x)QW\ao\,\ai\{z)) +{nx'^+eX+u) {u{x)Q}V\ao\,\ai\{z)) +{Q.X+eUj) {u{x)QW\ao\,\ai\{z)) =0 

(157) 

Simply, fll57p is 

u{x) ^ + + . + ^) = (158) 

From ( ]158p . we obtain 

2 TT dx + Cq ; Co, Ci = arbitrary constant (159) 

Q^\ao\,\a,\iz) 

Substitute (^M) into (^M) with Co= & Ci = 1. 

-iyp—{^fj.x^+ex) 

The (11601) is called " The first kind integral formalism of the second independent solution of 
the grand confluent hypergeometric function when u = negative integer, and 1." 
(ii) Case of u = positive integer except 1 

When u is greater than 1, A2 = l — u is negative integer. So A2 is smaller than Ai, by using 
the Case E) in theorem 1, gi{x) = G{x)\x=x^=i^u is equal to 'RyV'4,o,i(>i(^7 = |(1 + ^)',z = 
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^/ix^ ). Now, let assume that the second independent solution is 



g2{x) = u{x)gi{x) = u{x)TZW^q^^^ ( 7 = ^(1 + i/); 2; = -^/ix^ 



(161) 



By using similar process as we did before, 

g2{x) = 7^>Vv,o,Vl(^) 



dx 



(162) 



The fll62p is called "The second kind integral formalism of the second independent solution 
of the grand confluent hypergeometric function as = positive integer except 1." f ll60p and 
(11621) are inconvenient for direct calculation, as the integrand contains [QVV|ao|jQ,i|(z)] ^ 
and [T^W^o,^^ (z)] ^. Now let's investigate f ll60p and f ll62p more precisely and see how 
these two analytic functions have power series. These two functions can be described by the 
Frobenius method. 

(A)the case of u =0 or negative integer 

As is or negative integer, A2 = 1 — z^ is a positive integer. So Ai = is smaller than 
A2. By using Case C) in theorem 1, the first independent solution is 



gi{x) = G{x) 



A=Ai=0 



(163) 



The second independent solution is 

00 

g2ix) = gi{x)lnx + x^^ ^bn{\2] 



X 



(164) 



n=0 



To generate this solution, keep the recurrence formula in term of A and use it to find to 
coefficients Cn{n > 1) in terms of both A and Co, where coefficient Co remains arbitrary. 
Substitute these C„ into G{x) = x^ J2'^=o ^nx"' to obtain a function G(A, x) which depends 
on the variables A & x. From the Case C) in theorem 1, the second independent solution is 



92{x) 



dC{X,x) 



dX 



A=A2 



Inx ^ C„x" 



+A 



n=0 



+ J2 



A=A2 n=0 



dX 



(165) 



A=A2 



Plug C{x) = g{x)domin. + g{x) small into (|165|). 

dg{x) 

domin 



92[X) 



dX 



dg{x] 



small 



A=A2 



dX 



(166) 



A=A2 
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The first term of r.li.s in fll66p is 

dfjix) domin 



dX 



\=l-iy n=0 

oo 

lnx^C2„(A,n) 



njx 



2n+\ 



X=l-u 



X 



2n+X 



n=0 



X=l-u n=0 



gC2n(A,n) ^2n+A 



(167a) 



X=l-u 



Tfie second term of r.h.s in fll66p is 

dg{x) small 



dX 





I ^ C2n+l{X, 

n=0 

oo 

lna;^C2„+i(A,n) 



x=i-u 9X 



n)x 



2n+l+A 



A=l-i/ 



2n+l+A 



n=0 



A=l-i/ „=o 



^C'2n+l(A,n) 2n+l+A 

dX 



(167b) 

X=l-h' 



From f ITTll . 



Co — Co 



n-l 



C2n = CoY\ B2k+i ; > 1 



(168) 



fc=0 



Substitute i?„ in (!35|) into (I168p . and plug it into (1 16 Tap . 



dg{x] 



domin. 



dX 



X=l-v 



{l + \Yi{x))z'-'<A^,{r.z) 
(-1)-^,! (^0 + 1-7)! 



1 (h - . 



^ 2 n!(n + 1 - 7)!(^o - ^)! I (A: - ^o) (A: + 2 - 7) (A; + 1) 



n+l— 7 



(169) 



According to fl62l) . we have 



C3 = Co(Ao52 + A25i 

n— 1 



J- n— 1 n— 1 ^ n— 1 J— 1 \ n— 1 

C2n+1 = C'o^^O JJ^ -B2P+2 + ^ ( JJ^ -B2P+2 W B2k+1 j + ^2n Y\_ -^2p+l p > 2 

p=0 j=l ^ p=i fc=0 ^ p=0 ^ 

(170) 



From fl35l) . we get 

-Br) = 



£:(n + A + w) 



e(n + X + u) 



{n + l + X){n + X + u) (n + l + A)(n-l + 27 + A) 



fi{n -1 + X) + Q 



^{n -1 + X) + VL 



(n+l + A)(n + A + z/) (n + 1 + A)(n - 1 + 27 + A) 



(171a) 



(171b) 
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According to fll49ap and fll49b|) . put =-2^(^)0 + | - f ) = -2/i(V'o + 1 - 7) as n is odd 
in ffTTTb]) . As n is even, = -2/i(^^/'i + 1 - |j=-2/i(^^/'i + | - 7). Substitute f lTTU]) with 
f irTTaj) and fllTlbl) into fll67bD . then we obtain 

^72(x) = 7^WJ„,^, (7 = ^(1 + iy);z = —fix' 



z'-U {l + \n{x))A^,{r,z) 



-l)"l('o! (¥'0 + 1-7)! 



n-1 

E 



1 



. 2 n!(n + 1 - ^)\{^o -n)\^^{{k- ^0) (A: + 2 - 7) {k + 1) 



e r(^o + 2-7) (-^oW 

2"^ r(2-7) V^^Z^ 



n=0 



n!(2-7). 



ip\—n 



+ 1 - 7 + f )r(n + \)V{n + I - 7)(r2 - ^ 

X ^ . o , ^ : 2^ 



E 

fc=C 

+E 



r(A; + n + f )r(A; + n + f - 7) 



fc=0 

{-^o)nZ- A (r^ + l-i)(n + l-7 + |)^t^"^(r^ + i)^(n+|-7)(n 



n=0 

CXD 



n! (2-7)„ V2 



E^^n 



n=l 



(n + i)(n+|-7) 



n— 1 / , N 71—1 



E 

fc=0 



r(n + A; + |)r(n + A; + | - 7) 



E 



+ 



+ 1 - 7 + 



n— 1 / , N n— 1 



2(n + i)(n + |-7)^jL (p+l)(p + 2-7) 

00 1 / ■ , 1 I Lij \ n— 1 / , N 



iiJo -k) (A; + 2 - 7) (A; + 1) 



Npv- (j + l-7 + f) i-r 



X 



n=2 j: 



^2(j + |)(j + f-7) J-4(p + i)(p+i-7)fLV^ + ^)^^ + 2-7) 



i-1 

n 



(A; - 



(172) 



1 



+ 



+ 



(^0 - (/ + 2 - 7) (/ + 1) 



+ E 



1 1 

+ 7^ + 



1 



(^1-0 (^ + 1) (/ + I-7) 



(B)the case of = positive integer except 1 

As u is positive integer except 1, A2 = 1 — is negative integer. So A2 = 1 — is smaller 
than Ai = 0. By using Case E) in theorem 1, the first independent solution is 









1) 



h (7 = ^{^ + ^y,z = -^^^x') 

The second independent solution is 

00 

92ix) = ^i(x)lnx + x^' bniXi)x'^ 



(173) 



(174) 



n=0 



To generate this solution, keep the recurrence formula in term of A and use it to find to 
coefficients C„(n > 1) in terms of both A and Cq, where coefficient Cq remains arbitrary. 
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Substitute these C„ into G{x) = ^nx"" to obtain a function G(A, x) which depends 

on the variables A & x. From the Case E) in theorem 1, the second independent solution is 



92{x) 



dGi\,x) 



A=A 



Inx J2 Cna;"+^ + J2 



1 n=0 



da 



x=x 



1 n=0 



dX 



A=Ai 



(175) 



Plug G{x) = g{x)domin. + g{x) small into 01751). 



92{X) 



dg{x] 



domin. 



+ 



dg{x] 



small 



A=Ai 



A=Ai 



(176) 



The first term of r.h.s in fll76p is 



dg{x] 



domin. 



dX 



|-{$:C.„(A,n)x--^} 



A=0 - ■ n=0 

oo 



A=0 



lnx^C2„(A,n)a;2"+^ +^ 



ra=0 



dC2n{\n) 



A=0 n=0 



dX 



-X 



2n+A 



(177a) 



A=0 



The second term of r.h.s in fll76p is 



dg{x) 



small 



dX 



^{$^C^2„+i(A, 



n)x 



2n.+ l+A 



A=0 " - • - n=0 

oo 



A=0 



lnx^C2„+i(A, 



n)x 



2n+l+A 



n=0 



+5: 

A=0 n=0 



dC2n+liX,n) 2n+l+X 

dX 



(177b) 



A=0 



Substitute ffTTOj) . ffTTTal) and ffTTlbl) into ffTTTaj) and ffTTTb]) with A = as we did before. 
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Then, we obtain 



92{x) = QW(^,\,M {l = lil + ^y,z = —fix') 

n , w / N , rdaol +7) (-|ao|)n 
(1 + ln(x))i^„„| (7; z) + 



t 2 ^' (7)n 



fc=0 



(A;-|q;o|) + (A; + 7) 



— ,t;- 



g r(|ao|+7) 
■2 r(7) 

(-|«o|)n / 1 



l"o| 



(-|«o|)n J",^"(n + |)r(n+i)r(n + 7-|)(n-|ai|), , 



inx 



n=0 



n!(7). 



E 

A;=0 



+E 



(r! + l)(ii+|) 



|Qi|-n 



r(A; + n + |)r(A; + n + 7 + i) 
r(n + i)r(n + 7-i)(n-|ai|)fc , 



^n!(7)nV2 (ri+|)(n + 7-i)y ^ r(n + A; + |)r(n + A; + 7 + 1) 



2 



E TTfiTn 



(p- l«il) 



n-l 



n(7-i)i.V^^ + i)(^^ + ^ + i)fe^ V(l«i|-^) (^ + i)(^ + 7 + |) 



1 ^ (2(A; + l)+7) 



n + 



n-l 



(P- l"o|) 



rt— 1 

E 



2(n + + 7 - 5) (P + 1)(P + l) V (bol - *) (* + !)(<: + ?) 



1 + J^i±i±2L , 



00 n— 1 



^^it2(j + |)(j+7 

■i-i 



n-l 



irn 



(p- l"ll) 



i-1 

n 



(p+i)(P+5 + 7)ii(^ + l)(^ + 7) 



X 



(=0 



(|ao|-0 (Z + l)(Z + 7) 



n— 1 



1 ^ (2(A; + l)+7) 



(|ai|-0 (Z + |)(Z + i + 7) 



^"V (178) 



(C)the case of 1/ = 1 

As 1/ is 1, then two solution are identical: (Ai — A2). Then two solutions can be solved 

by Case B) in theorem 1, the first independent solution is 



gx{x) = G{x)\^^^^^^^^ Q>V|«o|,|ai|(7 = l]z = —^ix"^) 



(179) 



The second independent solution is 



92{x) = gi{x)\nx + x^' ^ 6n(Ai)a;" 



n=0 



(180) 
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And we can obtain its solution by using the similar process as we did before. 

1 



g2{x) = Q>Vf,„|,|,^| (7 = 1;^ = --/^x^) 



|ao| 



;-ir(i«oi!)^ 



n-l 



;2(n!)^(|«„|-n)!;^U*-|"ol) (* + 1) 



■4r(|.„i ^ .n.g izMv... i-i)^^ 



n=0 



fc=0 



(r(A; + n+|))2 



^^ {-\a,\U- (l (n + f) \ l"^" (r(n + |))^(n-|ai|), , 



n=0 

oo 



-E -n 



(n!)2 V- 



(^ + 5); fn^ (r(n + fc + |))2 



2 ^ / fe=0 
1 



+ 



l^2 11 rr,^ 1^2 



n-l 



2(n+|) 

oo n— 1 



+ 



J-1 

/=0 



1 2 

+ 



fc=0 
n— 1 



(i«oi-o (i+i); trV(i«ii-o (^+1) 



(A; + 1)2 
1 



(181) 



IV. INFINITE SERIES 



A. V — non-integer 



In the previous section, we discussed the solutions of polynomial case. The 1** indepen- 
dent solution has two eigennumber which are 



l«o| = -£^0,l,2,3, 



2n 

_ 1 
"2/i 2 



0,1,2,3, 



; only if |q;i| > |q;o| 



Also, the 2"^^ independent solution has two eigennumber which are 



(182) 



^o = -4+7-1^0,l,2,3,--- 

V'i = -S+7-|^0,l,2,3,--- ;onlyifV'i > V'o 
These conditions make the solutions as finite series. However, if 



(183) 



2n ' 



0,1,2,3,--- 



13 5 
2' 2' 2' 



2// / 



1-7,2-7,3-7,- 

2 ~ 7, 2 ~ 7, 2 ~ 7, 



(184) 
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Then, independent solutions turn to be infinite series. It means that as x goes to infinite, the 
function is divergent. Physically, this is means nothing. However, mathematically it can be 
an answer. From a mathematical point of view, the finite series is considered as the special 
case of infinite series. But it is totally a wrong approach. First of all, we need get a condition 
what situation makes polynomial as solution, after that, we need to expand it as infinite 
series. Polynomial case is totally different from the infinite series case. Polynomial case is not 
the special case of infinite series, they are independent with each other. Lets now consider 
the case of infinite series. From the 1'^* kind of grand confluent hypergeometric polynomial, 
let |ao| = and \ai\ = ^ | dSSD, f!69ap and fl69bl) where |ao|, |ai| 7^ 0, 1, 2, 3, ■ ■ ■ . 
Then the series turns to the infinite series for the case of Ai = 0. Also, from the 2"'^ kind of 
grand confluent hypergeometric polynomial, let i/jq = — — 1 + 'y and ipi = —-^ — ^ + 'y in 
f ll5ip - (ll53p where ipo, ipi 7^ 0, 1, 2, 3, ■ ■ ■ . Then the series turns to be the infinite series for 
the case of A2 = 1 — z^. Also, we skip these solutions because of too much details. 

B. V = integer 

The integral form of second independent solution for the case of z/ = negative integer, 0, 
and 1 is 



The fllSSp is called " The first kind integral formalism of the second independent solution of 
the grand confluent hypergeometric function as = negative integer, and 1 for the case 
of infinite series." Also, the integral form of second independent solution for the case oi u = 
positive integer except 1 is 



The fll86p is called "The second kind integral formalism of the second independent solution 
of the grand confluent hypergeometric function as v = positive integer except 1 for the case 
of infinite series." fllSSp and (11860 are inconvenient for direct calculation, as the integrand 



how these two analytic functions have power series. These two functions can be described 
by the Frobenius method. 

(A)the case of u =0 or negative integer 




(185) 




(186) 



contains [QW(2)] ^ and [7^>V(;^)] I Now let's 



investigate (11850 and (11860 more precisely 
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Replace i/jq and ■j/'i by — |^ — 1 + 7 and ~ | + 7 in fll72p . and change the maximum 
index number n from ipo to 00 in it. Also, change maximum index number k from ipi — n 
to 00. Then we obtain the second independent solution of grand confluent hypergeometric 
function as z/ = or negative integer for the case of infinite series. 

(B) the case of z/ = positive integer except 1 

Replace |ao| and \ai\ by and — — | in (11781) . and change the maximum index 
number n from |q;o| to 00 in it. Also, change maximum index number k from |a;i| — n to 00. 
Then we get the second independent solution of grand confluent hypergeometric function as 
u = positive integer except 1 for the case of infinite series. 

(C) the case of z/ = 1 

Replace |q;o| and \ai\ by and ~ | in (I18ip . and Change the maximum index 
number n from |q;o| to 00 in it. Also, change maximum index number k from — n to 00. 
And we obtain the second independent solution of grand confluent hypergeometric function 
as u = 1 for the case of infinite series. 

Furthermore, we will publish soon for the first time how to solve mathematical equations 
having three term recursion relations and go on producing the exact solutions of some of 
the well known special function theories that include Mathieu, Heun, Lame and the Grand 
Confluent Hypergeometric Function. We hope these new functions and their solutions will 
produce remarkable new range of applications not only in supersymmetric field theories 
as is shown here, but in the areas of all different classes of mathematical physics, applied 
mathematics and in engineering applications. 
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